Buffering plasmons in nanoparticle waveguides at the virtual-localized transition. 
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We study the plasmonic energy transfer from a locally excited nanoparticle (LE-NP) to a linear 
array of small NPs and we obtain the parametric dependence of the response function. An analytical 
expression allows us to distinguish the extended resonant states and the localized ones, as well as 
an elusive regime of virtual states. This last appears when the resonance width collapses and before 
it becomes a localized state. Contrary to common wisdom, the highest excitation transfer does not 
occur when the system has a well defined extended resonant state but just at the virtual-localized 
transition, where the main plasmonic modes have eigenfrequencies at the passband edge. The slow 
group velocity at this critical frequency enables the excitation buffering and hence favors a strong 
signal inside the chain. A similar situation should appear in many other physical systems. The 
extreme sensitivity of this transition to the waveguide and LE-NP parameters provides new tools 
for plasmonics. 
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I. INTRODUCTION 



Electromagnetic energy can be focused and guided be- 
low the light diffraction limit by transforming it into a 
collective plasmonic excitation, the surface plasmon po- 
lariton, that can propagate along a one dimensional array 
of nanoparticles (NPs) 1]. This feature has attracted sig- 
nificant attention due to its potential applications in op- 
toelectronic devices, sub- wavelength waveguides, random 
lasers, optical traps, and hot-spot based plasmonic sen- 
sors for ultrasensitive spectroscopy [J-Q . Previous works 
have already addressed the question of how to achieve a 
high degree of localization of plasmonic excitations Q 
and have studied the plasmon propagation on the waveg- 
uide formed by an ordered NP chain 0, 0]. However, a 
fundamental question remains open: how to transform a 
localized excitation into a strong signal somewhere else 
inside a finite chain. A natural idea, would be to exploit 
the divergent local density of states of 1-D systems, i.e. 
the vanishing group velocity at the band edge, as pro- 
posed for light buffering in photonic waveguides 6] . How- 
ever, such divergences are a bulk property absent in finite 
and semi-infinite chains, where the local density of states 
near the extremes cancels out at the passband edge. In 
a quantum tight-binding model, this corresponds to a 
semi-elliptical density of states. Therefore, the only al- 
ternative would seem to tailor the surface inhomogeneity 
to generate ad-hoc resonances Q that could easily be ex- 
cited and transfer energy. This poses a big challenge for 
both experiments and numerical simulations, since there 
is a broad range of parameters to be explored. Thus, this 
work resorts to a model that, containing only essentials, 
could be solved analytically by using a response function 
formalism. Besides the expected resonant and localized 
eigenmodes, we find an elusive regime of virtual states. 
This appears when the resonance width collapses and 



before becoming a localized state, providing for a conti- 
nuity between them. Quite surprisingly, we prove that 
virtual to localized states transition provides the route 
to optimal excitation transfer by recovering, at the chain 
extremes, a divergent local density of states with slow 
group velocity. The work is organized as follow: In Sec- 
tion II we present the model used to describe the system 
and develop a response function formalism and a pole 
analysis for it. Then, in Section HI we show and discuss 
the results, and finally, in Section IV we summarize the 
main conclusions of the work. 
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FIG. 1: (Color online) - An external source (Jilo) excites the 
surface plasmon (t^spo) of a NP which is coupled, through 
a^'^uj'i, to a NP waveguide, of bandwidth \lo 
where detection takes place. 
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II. LINEAR RESPONSE IN NANOPARTICLE 
WAVEGUIDES 

The system studied, see Fig. [l] is a linear array of 
N metal NPs coupled to a locally excited (LE) NP, de- 
scribed in the coupled dipole approximation (3, H, H, @| ■ 



The induced dipole moment Pi at r NP satisfies 
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Here, r^, wpi, wspi and rji correspond to the radius, bulk 
and surface plasmon frequencies, and electronic damping, 
eo is the free space permittivity, are 
respectively the external field and the electric field at 
the z*'' site produced by the j"* NP. In general, 1^ j,i 
is a complex function that depends on the separation 
vector ~(i j^i = dj idj i between NPs and the wave vector 

k. However, if d is small, ~^ j^i can be evaluated in the 
near field approximation: 
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where n is the refractive index of the host material. For 
a linear array of NPs, plasmon oscillations can only be 
transverse (T) or longitudinal (L) to the chain axis, and 
due to the cubic dependence of E on d, it is a good 
approximation to neglect contributions beyond nearest 
neighbors 0. Arranging all l^i and ^[^'■'as vectors P 
and E, Eq. [1] reads: 



lE =xE, 



where R is a diagonal matrix with i?^, 
and M is a tridiagonal matrix, with 



(3) 

-(4/3)7rrfw^,eo 



(4) 



for the LE-NP, while for any of the equidistant identical 
NPs along the chain 



~2 



: wgp - ir]uj, (Vi 7^ 0). 



The dipole-dipole coupling strengths are 
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with 7^ = 1, and 7^ — —2. Since we are interested in 
the particular case where only the LE-NP is different 
from the rest, it is convenient adopt the nearest- neighbor 
coupling as j — '^x ^^r i and j ^ 0. 

This description is accurate for: 1) fed ^ 1 and 2) 
r/d ^ 1/3, where higher order multipoles are negligible 
[10|. These conditions require small r's and hence a neg- 
ligible radiation damping correction [8j. 

Clearly, x is a response function (RF), hence excita- 
tion dynamics between the different sites i and j is de- 
termined by the corresponding matrix elements of x- In 
this way, the square dipole moment of the m^^ NP, 
produced when the LE-NP is externally excited with an 
electric field Eg, is: 
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Notice that (Iw^ — M) can be identified with a Green's 
function [HI 111]- The precise correspondence is presented 
in the Appendix and serves to exploit the analytical tools 
developed in this context. 



III. POLE ANALYSIS OF THE RESPONSE 
FUNCTION 

Since M is tridiagonal, Xij ('^^) admits exact analytical 



expressions as continued fractions |12l . |13[ . For finite sys- 
tems, Xij i'^'^) h^'S a set of isolated poles at the eigenval- 
ues of M whose real and imaginary parts are respectively 
the eigenfrequencies and their damping. Extending M 
to an infinite case enabled us to find a close expression 
for the response of our system at an arbitrary position m 
given an excitation at position i = 0: 



(8) 



where 



"1,0 = w|io/w|, (9) 
accounts for the surface asymmetry, while 

C\i^) = H^x/'n) ^ K ± ik (10) 

is a generalized wave vector. The term Xoo is the RF of 
the LE-NP, 
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where the factor 
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accounts for the effective coupling strengths, i.e. a = 
describes an isolated LE-NP. For small a, the peak at 
cjgpQ is further shifted and broadened by an(a;), a com- 
plex "self energy" accounting for the linear array. When 
N ^ 00: 



H(c.) = i[c.2_S 
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In the weak damping limit (WDL), i.e. 77 — > 0^, the 
propagating frequencies are given by the usual dispersion 
relation [2], 



w^(fc) — LUgp — 2Lu^cos{kd), 



(14) 



with wave number k € [— tt/c?, tt/c?]. Within the "pass- 



band" 
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< 2a;x, each frequency component of 



the excitation propagates with group velocity 
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Components outside the passband decay exponentially 
along the chain within the localization length k~^. The 
inclusion of electronic damping rj adds a further decay 
and smears out the dispersion relation. The overall be- 
havior of Eq|5] is consistent with the numerical solutions 
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including full retardation effects under similar conditions 
i- 

The local density of plasnionic states (LDPS) at site 
i = is given by Im(xpp). In the WDL, LDPS quanti- 
fies the participation of this site on the different eigen- 
frequencies in a range dui around u. More generally, for 
finite 7], ajlni(xoo) is proportional to the power absorbed 
when site is irradiated with frequency uj. 

The case a = of Eq. [TT] exemplifies the general be- 
havior of finite systems, where poles of the RP (zeros 
in the denominator of Xoo ) determine the frequencies of 
maximum energy absorption. In this situation, dissipa- 
tion occurs due to the damping processes accounted by rj. 
In an infinite system, a new mechanism appears as Im(n) 
also describes the irreversible energy spread through the 
chain. In this case, Eq. [TT] has "poles" , ujpoie , which 
solve: 
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The solution of this equation is: 
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The 



analysis of the different types of "poles" or solutions of 
Eq. [T7] is simplify in the WDL. In this case, two types 
of non-physical poles appear: 1) The pole has a positive 
imaginary part. 2) The pole is real but it corresponds to 
a divergence of the original equation (Xoo) with a non- 
physical self-energy (Im(n) < 0). This second alternative 
can appear because, in order to obtain a closed expres- 
sion for ujpoie, the denominator of Xoo should be squared 
(Eq. [T6| . This makes physical and non-physical self- 
energies indistinguishable. Therefore, both solutions are 
present in Eq. [TT) The second type of non-physical poles 
are associated with virtual states, because even though 
they do not appear as resonances, they aff'ect the LDPS 
within the passband, in the same way as localized or 
"real" states. 

Fig. [3] shows the real part of the poles of the response 
function, Wpoie: as function of a^/^ for four different val- 
ues of wspo ■ The figure also compare the case of the an- 
alytical solution resulting from the WDL of Eq. [17] with 
its numerical evaluation for a case with realistic damp- 
ing. Numerical evaluation is required for 77 because 
/3, V and Wp^/e depend on uj in this case. Notice the ac- 
curacy of the WDL approximation. These poles provide 
for discrete localized levels (L), resonant levels (R) and 
virtual states (V) according to the parametric region. 

In the WDL, when a pole is complex with a negative 
imaginary part, its real part lies within the passband and 
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FIG. 2: (Color online) Real part of poles of x ^ function 
of a^^'^. Here, Wx/i^ip = 0.25. In a) three different regimes 
are shown: discrete localized levels (L^ and Lu), resonant level 
(Ru) and virtual states (Vd and V^). Subscript u and d stand 
for "up" and "down", the possible positions of Re(ajpo;g) rela- 
tive to aj|p. The numerical solutions for finite r\ are compared 
with the analytical WDL. In panel b) we explore the depen- 
dence of the poles on ojspo- The colored dots indicate the 
parameters used in the discussion of the LDPS (see Fig. 3). 



corresponds to the eigenfrequency of the LE-NP (Fig. 
[3]) , while its imaginary part roughly represents the decay 
rate. This is the case of a resonant state. When the pole 
is real, the usual situation is that the system has a lo- 
calized eigenmode whose eigenfrequency lies outside the 
passband (Fig. [3]). An excitation at this frequency will 
remain indefinitely within the localization length 1/k. 
These two situations would typically exhaust the anal- 
ysis. However, for quantum systems it has recently be- 
come clear that the transition between these two regimes, 
although covering a very narrow parametric range, has 
subtle and unique properties: there is a real "pole" which 
nevertheless does not correspond to an eigenstate of the 
system (see [lj,ll5|)- As such, one might not know what 
to expect. This situation worsens in a plasmonic case 
where such parametric region broadens. Even though, 



4 



they are not physical poles, in fact they do not solve 
[w^ - Wgpo] - aU = 0, they still affect the LDPS, and 
hence the RF, within the passband (Fig. [3|) in a simi- 
lar way than localized states. In both cases, the LDPS 
is modulated at the band edges by a divergent factor 
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FIG. 3: (Color online) LDPS as function of the excitation 
frequency for different a values. Here, ojIpo/'^sp = 1-1 s^nd 
cjx/'^sp = 0.25. Also, Tj = Tjo — except for the panel show- 
ing a localized state case, where r) — 10~^. Notice the identi- 
cal shape of the LDPS near the passband edge for both virtual 
and localized states, and the divergence at the virtual-localized 
transition. 

This relationship between the poles of RF and the 
LDPS of the LE-NP, can be explicitly written in the 
WDL: 



LDPS 



Rqo 



X N (cj) X cx L (cj) . 



where = w^/w^p, c is a normalization constant 

2^2 



(18) 
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V4F2(l-a)-(l~/?)2' 
L is a Lorentzian function (in a plot as function of a;^), 

f2 
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and N (w) is proportional to the LDPS at the surface site 
of a semi-infinite chain of identical NPs, 
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In the resonant state regime, the maximum and the 
width of L (ui) are given by the real and imaginary part 



and 
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In the virtual and localized state regimes, these expres- 
sions are still valid but L (w) is no longer a Lorentzian 
function as F^ is imaginary. In this case, it is better to 
write the LDPS at the LE-NP as: 



LDPS =^ xN{io) 



(24) 
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Notice the divergences that appear at the virtual- 
localized transition, i. e. when ujpoie reaches the passband 
edge. These divergences will strongly favor the excita- 
tion transfer. Fig H] illustrates the frequency dependent 
excitation transfer along the chain for two cases: a) a 
system with a resonant state and b) a system with both, 
a virtual and a localized state. 

As occurs for the oscillating-overdamped transition, 
and many of the spectral bifurcations omnipresent in 
non-Hermitian Hamiltonians [l6j , the critical parameters 
of Eq. [T7] occur when the argument under a square root 
vanishes and they identify the phase transition in the 
dynamical behavior [17]. The transition resonant-virtual 
occurs for lm{ujpoie) — 0: 



2/3 + -(3'^ -I 
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and the transition virtual-localized takes place when 

(1-/3) 



ujpoie G 3? and w^^j^ = 1 ± 2^: 



ac2(±) = 2 ± 



V 



(26) 



These expressions, obtained in the WDL, will result es- 
sential to asses the role of the poles in the excitation 
transfer. 



IV. OPTIMAL EXCITATION TRANSFER 

Fig [5] shows the maximum excitation transfer, |Pm|2, 
enabled by a variation of w at each system configura- 
tion. Superposed are the critical values resulting from 
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FIG. 5: (Color online) Color scale shows the maximum exci- 
tation transfer to the S**" NP (max |P5(a;)p in arbitrary units) 
for 1] = O.OlcJsp and rjo = 0.02ajsp- Continuous line delim- 
its the region of resonant states (Ea l25|) and dash lines mark 
the virtual-localized transitions (Eq I26|) . both evaluated in the 
WDL. 
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FIG. 4: (Color online) Square dipolar moment (in arbitrary 
units) as function of the NP position and excitation fre- 
quency. Here, r; — O.Olusp, Vo = 0.02ajsp, '^ipo/'^sp = I-^j 
and uj^/oj^p = 0.25. In Fig.a a^/^ = 0.63 and in Fig.b 
a^/^ = 1.55. L, V, and R stand for localized, virtual and 
resonant states respectively, while u and d stand for up and 
down, the position of the pole relative to cjsp. 



Eqs. [55] and [211 Consistently with the above discus- 
sion, an appreciable transfer occurs in the resonant state 
regime. However, the maximum appears at the transi- 
tion between virtual and localized states. Note that the 
optimal configuration for excitation transfer does not oc- 
cur for cjspo = '^SP and a = 1, where the LE-NP is 
indistinguishable from the others, as one might naively 
expect. Instead, it occurs for the highly asymmetric con- 
figuration where a virtual- localized transition appears. In 
this case, an excitation of a local mode that is coupled 
to collective excitations at the band edge, has compo- 
nents with very slow group velocity. In consequence, un- 
der a continuous irradiation, the excitation can build up. 
This dynamical interpretation is consistent with the new 
strategy developed in the context of photonic crystals. 
There, waveguides with slow group velocities are used as 
a way of buffering light 6:]. Alternatively, one sees that 
just at the virtual- localized transition, the spectroscopic 
Eq. [5] favorably combines its three factors: the response 
function on the excited nanoparticle Xoo ('^) ^ high 
intensity peak, this peak is inside the passband, and the 

1/2 

relative effective coupling a^ Q is strong. Notice that a 



high transfer efficiency could be achieved even for wgpo 
so different from ujgp that the system could never form 
a resonant state just by changing a. In analogy with 
the addatom in an Anderson-Newns model , a strong 
interaction with the substrate captures a state from the 
continuum spectrum to build a second localized state that 
would constitute the "antibonding" orbital of a dimer 
[ist . This occurs through a virtual- localized transition 
and thus leads to the optimal excitation transfer shown 
in Fig [5] Experimentally, these critical points could be 
achieved by properly tuning the distance, radius, shape 
and material of the NPs. Additionally, this configuration 
acts as a very narrow filter for the external frequency in 
resonance with the passband edge (See Fig. [4]-b).The 
control of this critical phenomenon opens up many pos- 
sibilities for applications. For example, the extreme sen- 
sibility of excitation transfer on cJq.i when the system is 
close to a critical transition, would enable a new form of 
plasmon ruler suitable for biological and chemical appli- 
cations [l^l . This occurs because a varies with d^^ and, 
depending on wspo, wsp and which give the frequency 
offset of Fig [5l the system response sweeps through dif- 
ferent regimes within a narrow interval of a. Similarly, 
as small changes on the refractive index modify dramat- 
ically the coupling and hence the passband, the ex- 
citation transfer will also be extremely sensitive to the 
dielectric environment in a system tuned with the vir- 
tual-localized transition. 
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V. CONCLUSIONS 

We have demonstrated that, contrary to common wis- 
dom, the highest excitation transfer does not occur for a 
system with a well defined resonance but when a virtual 
state is transformed into a localized collective plasmonic 
mode whose eigenfrequency is just at the passband edge. 
The slow group velocity of an excitation with this criti- 
cal frequency enables the excitation buffering and hence 
favors a strong signal inside the chain. The extreme sen- 
sitivity of this transition to the waveguide and LE-NP 
parameters would provide new tools for plasmonics. As 
the basic model is quite general, our conclusions are uni- 
versal in nature and apply to any of the broad class of 
systems that can be mapped to a linear array of damped 
oscillators j20| . 
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APPENDIX. EQUIVALENCE BETWEEN 
PLASMONIC AND QUANTUM MECHANICAL 
MAGNITUDES. 



The correspondence between magnitudes in different 
models is analogous to the known Bloch discussion of 
unimpeded electronic motion in crystalline metals in 
terms of a periodic 1-d array of coupled pendula as de- 
scribed in Section 5.2 of ref. [11] and other textbooks. 
This is summarized in the following Table. 
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